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1. Introduction 



Gravitational-wave (GW) science has entered a new era. Experimentally Q, sev- 
eral ground-based laser-interferometer GW detectors (10-1 kHz) have been built 
in the United States (LIGO) [1], Europe (VIRGO and GEO) [2, 3] and Japan 
(TAMA) [4], and are now taking data at design sensitivity. Advanced optical 
configurations capable of reaching sensitivities slightly above and even below 
the so-called standard-quantum-limit for a. free test-particle, have been designed 
for second [5] and third generation [6] GW detectors (^ 2011-2020). A laser in- 
terferometer space antenna (LISA) [7] (10~^-10~^ Hz) might fly within the next 
decade. Resonant-bar detectors (^ 1 kHz) [8] are improving more and more 
their sensitivity, broadening their frequency band. At much lower frequencies, 
~ 10~^^ Hz, future cosmic microwave background (CMB) probes might detect 
GWs by measuring the CMB polarization [9]. Millisecond pulsar timing can set 
interesting upper limits [10] in the frequency range 10^^-10^^ Hz. At such fre- 
quencies, the large number of millisecond pulsars which will be detectable with 
the square kilometer array [11], would provide an ensemble of clocks that can be 
used as multiple arms of a GW detector 

Theoretically, the last years have been characterized by numerous major ad- 
vances. For what concerns the most promising GW sources for ground-based and 
space-based detectors, notably, binary systems composed of neutron stars (NS) 
and black holes (BHs), our understanding of the two-body problem and the GW- 
generation problem has improved significantly. The best-developed analytic ap- 
proximation method in general relativity is undoubtably the post-Newtonian (PN) 
method. The errors and ambiguities that characterized the very early literature on 
the PN problem of motion (for a review see, e.g., Ref. [12]), have been overcome. 
Robust predictions are currently available through 3.5PN order [13] (u^/c^), if 
the compact objects do not carry spin, and 2.5PN order [15] (w^/c^) if they carry 
spin. Resummation of the PN expansion aimed at pushing analytic calculations 
until the final stage of evolution, including the transition inspiral-merger-ring- 
down, have also been proposed, for the conservative two-body dynamics [16] and 
the radiation-reaction effects [17]. Quite interestingly, the effective-field-theory 
approach, commonly used in particle physics, has been extended to gravity, no- 

' GW experiments started with the pioneering work of Joseph Weber at Maryland in the 60s 
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tably to the two-body problem of motion [18]. The recent dazzling breakthrough 
in numerical relativity [19], with different independent groups being able to suc- 
cessfully evolve a comparable-mass BH binary throughout inspiral, merger and 
ring-down and extract the GW signal, is allowing to dig out details of the nonlin- 
ear dynamics which could not be fully predicted with other means. 

Our knowledge has also progressed on the problem of motion of a point par- 
ticle in curved spacetime when the emission of GWs is taken into account (non- 
geodesic motion) [20, 21]. To solve this problem is of considerable importance 
for predicting very accurate waveforms emitted by extreme mass-ratio binaries, 
which are among the most promising sources for LISA [22]. 

The GW community working at the interface between the theory and the ex- 
periment has provided templates [16, 17,23] for binaries and developed robust 
algorithms [24, 25] for pulsars and other GW sources observable with ground- 
based and space-based interferometers. The joined work of data analysts and 
experimentalists has established astrophysically significant upper limits for sev- 
eral GW sources [26-28] and is now eagerly waiting for the first detection. 

These lectures were envisioned to be an introductory, basic course in GW the- 
ory. Many of the topics that we shall address are thoroughly discussed in several 
books [29-35] and proceedings or reviews [36^2]. The lectures focused more 
on binary systems, probably because biased towards' the author own background 
and expertise. The lectures are organized as follows. In Sec. |2]we start by de- 
riving the wave equation in linearized gravity and discuss the main properties 
of GWs. In Sec. [3] we describe the interaction of GWs with free test particles 
and the key ideas underlying the functioning of GW detectors. Section|4]reviews 
the effective stress-energy tensor of GWs. Section |5] is devoted to the genera- 
tion problem. We explicitly derive the gravitational field at leading order assum- 
ing slow-motion, weak-gravity and negligible self-gravity. We then discuss how 
those results can be extended to non-negligible self-gravity sources. As a first 
application, in Sec.|6]we compute the GW signal from binary systems. We dis- 
cuss briefly the state-of-the-art of PN calculations and NR results. As an example 
of data-analysis issues, we compute the GW templates in the stationary-phase- 
approximation (SPA). In Sec.|2]we apply the results of Sec.|6]to other astrophys- 
ical sources, notably pulsars and supernovae. Section [8] focus on cosmological 
sources at much higher red-shift 2^1. We review the main physical mecha- 
nisms that could have produced GWs in the early Universe, notably first-order 
phase transitions, cosmic and fundamental strings, and inflation. 
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2. Linearization of Einstein equations 

In 1916 Einstein realized the propagation effects at finite velocity in the gravi- 
tational equations and predicted the existence of wave-like solutions of the lin- 
earized vacuum field equations [43]. In this section we expand Einstein equations 
around the flat Minkowski metric derive the wave equation and put the solution 
in a simple form using an appropriate gauge. 

2.1. Einstein equations and gauge symmetry 
The Einstein action reads 

S. = Y^ld'^V-9R, (2.1) 

where c denotes the speed of light, G the Newton constant, g^i, is the four dimen- 
sional metric and g = dct(g^^). Henceforth, we use the following conventions. 
The flat Minkowski metric is 77^1, = (— , +, +, +), Greek indices denote space- 
time coordinates /i, = 0, 1, 2, 3, whereas Latin indices denote spacelike coor- 
dinates i,j = 1,2,3. Moreover, = (a;°,x) = (ct,x), thus d'^x = cdtd^x. 
Partial derivatives 9^ will be denoted with a comma, while covariant derivatives 
with a semicolon. The scalar tensor R in Eq. ( 12.1b is obtained from the curvature 
tensor as 

R = g^^ Rfii, , R^u = gf"^ Rp^av , (2.2) 
where T'^,„ are the affine connections 

per 

^p ^ 1 ma , ^ _ ^\ n 4^ 

2 ^ \dxP dx- dx^ ) ■ ^ ' 

The curvature tensor satisfies the following properties 

Rpupa Ruppa Rpyap ; Rpupa Rpapu ; (2.5) 

Rpupa + Rpavp + Rppau = j R pvp\(j ^" R pav^p ^" R ppcr;u ^ ^ ■ 



(2.6) 

The latter equation is known as the Bianchi identity. We define the matter energy- 
momentum tensor Tpi, from the variation of the matter action Sm under a change 
of the metric g^^ gpi, + Sgp^, that is 

SS^^^ I d^x^T^" 5gp,. (2.7) 
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The variation of the total action S — Sg + Sm with respect to g^i, gives the 
Einstein equations 



1 „ SttG 
2- 



Gfiu — R^ij — — g^i, R — — — Tpjy . (2.8) 



The above equations are nonlinear equations with well posed initial value struc- 
ture, i.e. they determine future values of g^j, from given initial values. Since 
/i = 0, • • • 3, = 0, • • • 3, Eq. ( 12.8b contains sixteen differential equations, which 
reduce to ten differential equations if the symmetry of the tensors G^i^ and T^i, 
is used. Finally, because of the Bianchi identity we have G^^'"' = 0, thus the ten 
differential equations reduce to six. 

General relativity is invariant under the group of all possible coordinate trans- 
formations 

x^-^a;'^(x), (2.9) 

where x'^ is invertible, differentiable and with a differentiable inverse. Under the 
above transformation, the metric transforms as 

a^..{x) ^ (j'^ui^') ^Q^Q^ 9pAx) ■ (2.10) 

We assume that there exists a reference frame in which, on a sufficiently large 
spacetime region, we can write 

gtiu^Vi^i' + htiu, |ft.^i.|<l. (2.11) 

By choosing this particular reference frame, we break the invariance of general 
relativity under coordinate transformations. However, a residual gauge symmetry 
remains. Let us consider the following coordinate transformation 

a;'^^x"^ = x^+e^(a;), 19^^.1 < | VI • (2-12) 

The metric transforms as 

9^ = Vt.u - a.^M - d^i, + h^, + 0{de) , (2.13) 

thus, introducing 

^fj,u — h^i, ^^,1/ ?i^,At I (2.14) 

we have 

In conclusion, the slowly varying coordinate transformations ( 12.121 ) are a sym- 
metry of the linearized theory. Under a finite, global (x-independent) Lorentz 
transformation 

x^-^A^x-", APAZr^p^^r,^,, (2.16) 



Gravitational waves 



9 



the metric transforms as 

Sm'^ = A;:A^5p, = ?7^, + A^A^ V(a^) = + (2-17) 

thus, hf^i, is a tensor under Lorentz transformations. It is straightforward to prove 
that hfj.1, is also invariant under translations. In conclusions, linearized theory is 
invariant under the Poincare group and under the transformation — > .t^ + 

with la^e^i < 1. 

2.2. Wave equation 



Let us now linearize Einstein equations posing = 77^^ + /i^^. At linear order 
2 ( 

1 



in hfj,i, the affine connections and curvature tensor read 



= ^'f^ (^p^Am + 9^.hxp - dxKp) , (2.18) 



pp 2 

fipcr ^p^ p,(j 

more explicitly 



Kpa = a^r;!, - d.K^ + oih') , (2.19) 



R^upa — 2 i^pi^^pty ^~ ^ap^vp ^pp^i^a ^av^p^p^ • (2.20) 

Using the above equations, it is straightforward to show that the linearized Rie- 
mann tensor is invariant under the transformation hp^ — > hp^, — dp^i, — d,y^p. 
Equation (I2.20l i can be greatly simplified if we introduce the so-called trace- 
reverse tensor 

V^h^'^-^ifh, (2.21) 

where h ~ rjap h"'^ and h = —h, which explains the name. Some algebra leads 
to 

nK,+ij,,dP d%,x-dP dJip^^dP djip. + Oih') = ^T,,^ (2.22) 

where the wave operator □ = ripcd''d'^. To further simplify Eq. (12.221 ) we 
can impose the Lorenz gauge (also denoted in the literature as harmonic or De 
Donder gauge) 

d^h^'' = , (2.23) 

and obtain 

aK, = -^^n,. (2.24) 

If /i^'' does not satisfy the Lorenz gauge, i.e. S^/i'" = Qi,, we can introduce a 
coordinate transformation such that h'p^ = hp^ — ^p,^ — ^^.^ + rip^{dp^'') and 
impose D^^, = q^^, obtaining dph'^'^ = 0. 
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Summarizing, the Lorenz gauge imposes 4 conditions that allow to reduce the 
10 independent components of the 4x4 symmetric tensor /i^^ to 6 independent 
components. Note that we also have the condition d^T^'^ = 0, which is the 
conservation of the energy-momentum tensor of the matter in linearized theory. 
By contrast in the full theory T^^^ = 0. 

2.3. Tmnsverse-tmceless gauge 

We want to study the propagation of GWs once they have been generated. We 
set T^i, = in Eq. ( 12.24b and obtain the wave equation in vacuum 

Uhf,^ = . (2.25) 

GWs propagate at the speed of light. Within the Lorenz gauge we can always 
consider coordinate transformations such that ~ 0. The trace-reverse tensor 
transforms as /i'^^ = h^^ + with = rj^y dp^P - ^^^^ - Using 
^S.iiv = 0, we can subtract 4 of the 6 components of h^i,. More specifically, we 
can choose such that h ~ Q and ^' such that /t,*" = 0, thus OqIi'^'^ ~ 0. The 
GW being a time-dependent field, we can set /i"" = 0. We denote the field htj 
which satisfies the following transverse and traceless gauge conditions, 

h"^ ^0, /i"* = , d,h'^ = , h" ^ , (2.26) 

the transverse-traceless tensor hjl^ . Note that for a single plane wave with wave 
vector k and propagation direction n = k/fc, the transversality condition reduces 
to n' hj^ = 0. Without loosing in generality, we can assume that the plane wave 
propagates along the z-axis, thus 

I h+ \ 

hj^{t,z)~ hx -/^+ cos 
\ / 

where we indicate with and hx the two independent polarization states. Fol- 
lowing [31, 35], we can introduce the projector operator Pij{n) ~ Sij — Ui nj, 
which satisfies the conditions 

= Pj, , n'' = , P''' = P.t , P,;. = 2 , (2.28) 

and the A-operator 

A,, ki (n) - P,fc Pji-^ P^J Pki , (2.29) 
and obtain the TT field for a generic propagation direction 

hj^ = A,,- , (2.30) 



LJ [t 



(2.27) 
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where hki is given in the Lorenz gauge but not necessarily in the TT gauge. 

The GW is described in the TT gauge by a 2 x 2 matrix in the plane orthogonal 
to the direction of propagation n. If we perform a rotation t/j about the axis n, we 
obtain 

h^±ih+~^e'^^''^ihy,±ih+). (2.31) 

In particle physics we call helicity the projection of the (total) angular momentum 
along the propagation direction: 7i = J • n = S • n, being S the particle's spin. 
Under a rotation about the propagation direction the helicity states transform 
as /i — > e*'^ ^ h. Thus, Eq. ( 12.31b states that — ih+ are the helicity states and 
that the graviton is a spin-2 particle. 



3. Interaction of gravitational waves with point particles 

3.1. Newtonian and relativistic description of tidal gravity 

Let us consider two point particles, labeled A and B, falling freely through 3- 
D Euclidean space under the action of an external Newtonian potential $. We 
assume that at time t = the particles are separated by a small distance ^ and 
have initially the same speed 'vj^{t = 0) = vb(< = 0). Since the two particles 
are at slightly different positions, they experience a slightly different gravitational 
potential <I> and a different acceleration g = — V<i>. At time t > {), VA{t) ^ 
VB(t). Let us introduce the separation vector ^ = x^i — in 3-D Euclidean 
space. We have 



rf2|^_ / \ 



e~£]i\ (3.1) 



where the second equality is obtained by Taylor expanding around the position 
of particle A. The tensor f j is called the Newtonian tidal- gravity tensor [42], it 
measures the inhomogeneities of Newtonian gravity. It is the tensor responsible 
of the Moon's tides on the Earth's ocean. 

We now generalize the above Newtonian discussion to Einstein theory. In 
general relativity, nonspinning test particles move along geodesies 

d^x^' ^„ , , dxP dx"^ „ 

where r^g.(.T) is given by Eq. ( |2.19t . Let us consider two nearby geodesies, 
labeled A and B, parametrized by x^{t) and x'^{t) + ^^(t), with |^'^| smaller 
than the typical scale on which the gravitational field varies. By expanding the 
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geodesic equation of particle B around the position of particle A, and subtracting 
it from the geodesic equation of particle A, we get 

V.V.e^-«'V«'^^, (3.3, 

= dx^ /dr being the four- velocity and where we introduce the covariant 
derivative along the curve x^^ (r) 

V-e^-^^ + r^^.C"^- (3.4) 

Thus, two nearby time-like geodesies experience a tidal gravitational /orce pro- 
portional to the Riemann tensor 

3.2. Description in the transverse-tmceless gauge 

In this section we describe the interaction of a GW with a point particle in the TT 
gauge. Let us consider a test particle A at rest at time r = 0. Using the geodesic 
equation, we have 

dV _ _( , dx" dX^\ _ _ / ; dx'' dX^\ 

dr^W^o^ [ dT dr l ~ I ° rfr dr l ' ^^'^^ 

^ ' \t = Q ^ ' |t = 

Because the particle is initially at rest {dx^^ /dT)T-=Q = (c, 0) and 

Too - I V'' {dohoj + dohjo - d,h„o) ■ (3.6) 

In the TT gauge hoQ — and hoj = 0, so (rQQ)^=o ~ 0. Thus, we conclude that 
in the TT gauge, if at time r = 0, dx^ /dr = 0, also (P'x^ jdr^ = and a particle 
at rest before the GW arrives, remains at rest. The coordinates in the TT gauge 
stretch themselves when the GW arrives so that the coordinate position of the 
point particles, initially at rest, does not vary. What varies is the proper distance 
between the two particles and physical effects are monitored by proper distances. 
For a wave propagating along the z-axis, the metric is [see Eq. ( 12. 271 )1 



ds^ = -c^ dt^ + d'? + dy^ 



2; 

1 — COS uj \i 



1 + /i+ cos UJ \t 



z 
c 

If particles A and B set down initially along the a;-axis, we have 



2dxdyhx coslu (^t ^ .(3.7) 



dx' 



2 

where L is the initial, unperturbed distance between particles A and B 



L [1 + ^ cosujt I , (3.8) 
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3.3. Description in the free-falling frame 

It is always possible to perform a change of coordinates such that at a given 
space-time point Q, we can set T'^^{Q) ~ and {dP x^^ / (It"^) q = 0. In this 
frame, at one moment in space and one moment in time, the point particle is free 
falling (FF). This frame can be explicitly constructed using Riemann normal co- 
ordinates [30]. Actually, it is possible to build a frame such that the point particle 
is free-falling all along the geodesies using Fermi normal coordinates [30]. 

Let us introduce a FF frame attached to particle A with spatial origin at x^ = 
and coordinate time equal to proper time = t. By definition of a FF frame, 
the metric reduces to Minkowski metric at the origin and all its derivatives vanish 
at the origin, that is 



ds^ = -c^ dt^ + (ix^ + O [j^j ' (3-9) 

where 7?. is the curvature radius 7?.^^ = \Rfj,upa\- Doing explicitly the calculation 
at second order in x, one finds [44] 

ds^ = -c^ dt^ [l + RiOjQ x' x^] - 2c dt dx' ( ^-Rojife x^ 



dx^ dx^ 



Sij - -Rijki x'' x^ 



(3.10) 



For GW experiments located on the Earth, the interferometer is not in free fall 
with respect to the Earth gravity. The detector is subjected to an acceleration 
a = — g with respect to a local inertial frame and it rotates with respect to local 
gyroscopes. Thus, in general the effect of GWs on point particles compete with 
other effects. We shall restrict our discussion to the frequency band (10-10'^ Hz) 
in which the other effects are subdominant and/or static. 

Let us compute the equation of geodesic deviation in the FF frame attached to 
particle A. We have 

where in the last equality we use F"^ = 0. Since we assume that the particles are 
initially at rest, ~ Sq. Using = and the fact that F^j. q can be neglected 
when computed at position A, we have 



V„V„e = -^ ^ -^^-R^oko^. (3.12) 



To complete the calculation we need to evaluate the Riemann tensor R^qi.q- As 
discussed in Sec. 12.21 in linearized theory the Riemann tensor is invariant under 
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Fig. 1 . We show how point particles along a ring move as a result of the interaction with a GW 
propagating in the direction peipendicular to the plane of the ring. The left panel refers to a wave 
with + polarization, the right panel with X polarization. 



change of coordinates, so we can compute it in the TT gauge. Using Eq. (|2.20| i. 
we obtain 

-^Jofco = ^2c2 ■ (3.13) 



Thus, 



In conclusion, in the FF frame the effect of a GW on a point particle of mass m 
can be described in terms of a Newtonian force Fi = {iti/2) hjj^ . Note that in 
the FF frame, coordinate distances and proper distances coincide, and we recover 
immediately Eq. ( 13.8b . 

The description in the FF frame is useful and simple as long as we can write 
the metric as g^^, = rj^^, + 0(.t^/7?.^), i.e. as long as we can disregard the 
corrections /'R?. Since 7?."^ ~ \Rmjo\ ^ h ^ h/X^^, we have x^/TZ"^ ~ 
L^/i/Aq^, and comparing it with 5L/L ^ h, we find Li^/X^^ <C 1. This 
condition is satisfied by ground-based detectors because L ^ A km and Aqw ~ 
3000 km, but not by space-based detectors which have i ~ 5 x 10^ km and will 
observe GWs with wavelength shorter than L. [For a recent thorough analysis 
and a proof of the equivalence between the TT and FF description, see, e.g., 
Ref. [45].] 



3.4. Key ideas underlying gravitational-wave detectors 

To illustrate the effect of GWs on FF particles, we consider a ring of point parti- 
cles initially at rest with respect to a FF frame attached to the center of the ring, 
as shown in Fig.[T] We determine the motion of the particles considering the + 
and X polarizations separately. If only the + polarization is present, we have 

hj^ = h+ { \ \ ) sinwt , = [xq + 5x{t), yo + 6y{t)] , (3.15) 
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where Xq and j/q are the unperturbed position at time t = 0. Thus 



Sx{t) 



xq sinojt Sy{t) 



■ hq sin Lut . 



(3.16) 



If only the x polarization is present, a straightforward calculation gives 



Sx(t) = j/o sinijji 



5y{t) = Xq sintjt . 



(3.17) 



The + and x polarizations differ by a rotation of 45°. In Fig.|2]we show the lines 
of force associated to the + and x polarizations. 

The simplest GW detector we can imagine is a body of mass m at a distance L 
from a fiducial laboratory point, connected to it by a spring of resonant frequency 
fi and quality factor Q. Einstein equation of geodesic deviation predicts that the 
infinitesimal displacement AL of the mass along the line of separation from its 
equilibrium position satisfies the equation [36] (valid for wavelengths ^ L and 
in the FF frame of the observer at the fiducial laboratory point) 



AL(t) + 2^ AL(t) + AL{t) = - 



F+h+it) + F^h,,{t) 



(3.18) 



where -F+.x are coefficients of order unity which depend on the direction of the 
source [see Eqs. ( |6.31t , ( |6.31| l below] and the GW polarization angle. 

Laser-interferometer GW detectors are composed of two perpendicular km- 
scale arm cavities with two test-mass mirrors hung by wires at the end of each 
cavity. The tiny displacements AL of the mirrors induced by a passing GW 
are monitored with very high accuracy by measuring the relative optical phase 
between the light paths in each interferometer arm. The mirrors are pendula with 
quality factor Q quite high and resonant frequency O much lower (~ 1 Hz) than 
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f(Hz) 

Fig. 3. We plot the square root of the noise spectral density versus frequency for the three LIGO 
detectors together with the LIGO noise curve at designed sensitivity . The noise curves refer to June 
2006, during the fifth scientific run. 



the typical GW frequency {^-^ 100 Hz). In this case Eq. ( 13.18b . written in Fourier 
domain, reduces to AL/L ~ h. The typical amplitude, at 100 Hz, of GWs 
emitted by binary systems in the VIRGO cluster of galaxies 20 Mpc distant), 
which is the largest distance the first-generation ground-based interferometers 
can probe, is ~ 10^^^. This means AL ^ 10^^^ m, a very tiny number It 
may appear rather discouraging, especially if we think to monitor the test-mass 
motion with light of wavelength nearly 10^^ times larger. However, this precision 
is currently be demonstrated experimentally. 

The electromagnetic signal leaking out the interferometer's dark-port contains 
the GW signal but also noise — for example the thermal noise from the suspen- 
sion system and the mirror itself, can shake the mirror mimicking the effect of 
a GW. The root-mean-square of the noise is generally expressed in terms of the 
noise power per unit frequency Sh through the relation h ^ \/ Sh{f) Af ^ 
AL/L, AL being the mirror displacement induced by noise and Af the fre- 
quency bandwidth. In Fig.[3]we plot the noise curves of LIGOs (June 2006). The 
interferometers are currently operating at design sensitivity for almost the entire 
frequency band. 
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4. Effective stress-energy tensor of gravitational waves 

Until now we have defined the GWs as fluctuations of a flat spacetime. Here, 
we want to be more general and consider GWs as perturbations of a generic 
background g^^, that is 

9tJ.iy =9^,y + , \hf,^\<^l. (4.1) 

We need a criterion to define what is the background and what is the perturbation. 
Heenceforth, we follow closely the derivation in Ref. [35]. In general there are 
two cases: 

1. has typical scale Lb and h^^, has typical wavelength A with A ^ Lb, i.e. 
h^,y is a small ripple on a smooth background; 

2. has frequencies only up to Jb and /i^^ is different from zero around / 
with / ^ /b, i.e. the background is static. 

Let us expand through 0{h?). Note that we have now two small parameters 
h and A/Lb (or /b / /), we have 

= 2^ + R^^} + 4'-^ + • • • • (4.2) 

low frcq high frcq low and high frcq. 

Quantities having a bar are computed using the background metric 'g^^', they 
contain only low-frequency modes. The superscript (1) [(2)] in Eq. ( 14.21 ) refers to 
quantities computed at linear (quadratic) order in h. Using the Einstein equations 
we get 



low frcq 



(4.3) 



We introduce a scale I such that A <C ^ <C Lb, and average over a spatial volume 
[35,46]. We denote the average as (). Short-wave modes average to zero, 
whereas modes with wavelength Ls are constant. We can rewrite Eq. ( 14.31 ) as 

R,. = -«'> + ^(r,.-ig,.r) 

- -(4'.') + ^ . (4.4) 

Defining the effective stress-energy tensor of GWs 

V = {Ri'^ - Iv^u R^'^) , (4.5) 
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we have 

1 SttG 
R^u-^9^,.R=—{T^u + t^,). (4.6) 

An explicit calculation carried on far from the source gives 

(A 

t^^-^ = {d^Kp d^h'^f') . (4.7) 

For a plane wave, using the TT gauge 



^oo^^(A-/.->.^(.^ + .^>, (4.8) 



and the GW energy flux per unit area is 

dE 



dtdA IGnG 



{ill + hi). (4-9) 



For a supernovae dE/{dtdA) ^ f'^h? / {IQ-kG) ^ 400erg/(cm^sec), where 
we set /i = 10^^^ and / = 1 kHz. The GW burst has a duration of a few msec. It 
is telling to compare it with the neutrino energy flux ~ lO^erg/ (cm^sec) and the 
photon energy flux (optical radiation) ^ 10^"'crg/ (cm^sec) from a supernovae. 
Neutrinos and optical radiation are emitted during a few seconds and one week, 
respectively. 



5. Generation of gravitational waves 

5.1. Sources in slow motion, weak-field and negligible self-gravity 

In this section we evaluate the leading order contribution to the metric pertur- 
bations under the assumption that the internal motions of the source are slow 
compared to the speed of light. We also assume that the source's self-gravity is 
negligible. Henceforth, we shall discuss how to extend those results to sources 
with non-negligible self-gravity. We start from 

□V = -^^T^M-- dtX" = 0, 9^T^'' = 0, (5.1) 
and introduce retarded Green functions 
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which satisfy O^G{x — x' X — x'). The solution of Eq. ( 15.11 ) can be 

written as 

V(2;) = -^^ j d^x'G{x-x')T^,{x'). (5.3) 
Outside the source, using the TT gauge, we have 

7^J^(t,x) = A,,, fc^n) ^ j d^x' (^t - ^^^;x') . (5.4) 

Denoting by d the typical size of the source, assuming to be far from the source, 
i.e. r ^ d,we can write |x — x'| = r — x' • n + 0{d^ /r), and Eq. ( 15.3b becomes 



Tly. (i, x) ~ - — A,,- (n) / d^^' Tki[t-- + ; x' ) . (5.5) 



r c 



We can simplify the above equations if we assume that typical velocities inside 
the sources are much smaller than the speed of light c. If w is the typical fre- 
quency associated to the source motion, typical source velocities are v ^ uj d. 
As we shall see in the following, the GW signal is determined by the leading 
multipole moments, thus ijJgw ^ ^ v / d and Aqw ~ (c/^') d- If w/c ^ 1, we 
have Agw > d. 

Applying a Fourier decomposition, we can write 

(5.6) 

using ojx' • n d/c <C 1, expanding the exponential and Taylor-expanding 

Tki we get 



f d'xT'''{t,^)x"'+ 
c dt J 

^n„np^y d^xT'''{t,x)x^ xP + 



(5.7) 



The above expression is valid in linearized gravity and for negligible self-gravity 
sources, i.e. for sources whose dynamics is not determined by gravitational 
forces. We notice that in Eq. ( 15.7b the higher multipoles are suppressed by a 
factor v/c. To make Eq. ( 15.7b more transparent, we can express the momenta T^^ 
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in terms of the momenta of T"" and T°*. Let us first introduce the momenta of 
the mass density 

M = ^ J d^xT^"{t,x) , (5.8) 

M' = ^ J d^xT°°{t,x)x\ (5.9) 

W = ^ J d^xT°"{t,x)x' x\ (5.10) 

(5.11) 

and impose the conservation law df^iT^^ = valid in linearized gravity. Setting 
= 0, we have d^T^^ + diT"^^ = 0, integrating this equation in a volume 
containing the source, we obtain the conservation of the mass M = 0. Similarly, 
one can prove the conservation of the momentum M'- = Moreover, we have 



/ d^xx'x^doT"" = - [ d^xx'x^dkT^^ 
Jv Jv 

d^xix^T'^' + x'T^^), (5.12) 



where the second line is obtained after integrating by parts. Finally, 

M'^ =2 [ d^xT'^ . (5.13) 
Jv 

Thus, the leading term in Eq. ( I5.7l i. can be rewritten as 

hf;^it,^) = 1 ^ A,,,.,(n)M'=' (t- n , (5.14) 
r c V c/ 

where A/*^' is given by Eq. ( 15.111 ). The quantity T'^^ /c^ in Eq. ( 15.111 ) is a mass 
density. Besides the rest-mass contribution, it can contain terms due to the kinetic 
energy and the potential energy. For sources having strong gravitational field, as 
NSs and BHs, the mass density can depend also on the binding energy. Only for 
weak fields and small velocities, which is the assumption so far made, r"°/c^ 
reduces to the rest-mass density p. 

Henceforth, we shall discuss some applications to binary systems and pulsars, 
so it is convenient to compute explicitly /i+ and h^- Assuming that the GW 
propagates along the direction n = (cos0 sin^jSini/) sin cos a straight 
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calculations gives: 
G 



J < Mil (sin^ (f) - cos^ Q cos^ 0) + M22 (cos^ - cos^ 
-A/as sin^ Q - Mii sin 20 (1 + cos^ B) + M13 cos sin 26*+ 
M23 sin 26* sin^l , (5.15) 



2G f 1 •• 

— 7 i 7: (Mil - M22) cos 6* sin 20 - M12 cos0 cos 20 
re* 2 



Ml 3 sin0 sin0 + Mo-i cos sin 



in6l| . (5.16) 



5.2. Sources in slow motion and weak-field, but non-negligible self-gravity 

As already stated above, the derivation in linearized gravity of the quadrupole 
formula ( 15.141 ) cannot be applied to systems like binary stars whose dynamics 
is dominated by gravitational forces. In fact, because of the conservation-law 
valid in linearized gravity d^T^„ — 0, the two bodies move along geodesies in 
Minkowski spacetime. The extension to the case in which self-gravity is non- 
negligible can be done as follows [47]. 

In full general relativity one can define the field H^^ such that 

g^'' ^ - (5.17) 

where in the weak-field limit H^'^ coincides with the reverse-trace tensor that we 
introduced above. Imposing the harmonic gauge d^H^^'^ = 0, one derives 

= ^ [(-,9) T^,, + T^,] , (5.18) 

where r^i^ is the pseudotensor depending on i/^^ that can be read explicitly from 
Refs. [14,29]. The conservation law reads in this case 

5M[(-.9)r^" + Tn =0. (5.19) 



We can redo the derivation in linearized gravity (see Sec. |2.2t . but replace T^^ 
{—g) + T^i,, obtaining for the leading term in Eq. (|5.14| i J Tqo d^x — > 
/ (2oo + Too) a;' x^ d^x. For sources characterized by weak gravity tqo is negli- 
gible with respect to Tqo. Even though at the end one obtains the same formula, 
it is crucial to take into account the second order corrections in h^^, i.e. the field 
T^j^, in the conservation law. Otherwise the sources would be obliged to move 
along geodesies in Minkowski spacetime, instead of moving in a bounded orbit. 
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5.3. Radiated energy, angular momentum and linear momentum 

Using the results of Sec. |4] we can compute the power radiated at leading order 

„2 3 (-f 

where we introduce the traceless quadrupole tensor 

Q,j = -^S,jMkk- (5.21) 
Using the following relations 

dfl 1 , 

— Hi Uj = - Sij , (5.22) 

^ Tii Uj UkUi^ — {5ij 5ki + 6-ik Sji + 5ii Sjk) , (5.23) 
we derive for the total power radiated 

In the literature Eq. ( |5.24t is generally denoted as the quadrupole formula. 

GWs not only carry away from the source the energy, but also angular mo- 
mentum and linear momentum. At leading-order the angular-momentum radiated 
is [35] 

(5.25) 

while the hnear momentum radiated is given by [35] 
dP' G f . --TT 



dt Sttc' 



/■■ - TT 
du;QjkdVj;f. (5.26) 



Under parity, x — x, the mass quadrupole does not vary, and the integral in 
Eq. ( 15.26b is overall odd and vanishes. The first nonzero contribution comes at 
order 0{l/c^) from the interference between the mass quadrupole and the sum 
of the octupole and current quadrupole. As a consequence of the loss of linear 
momentum through GW emission, the BH formed by the coalescence of a BH 
binary can acquire a kick or recoil velocity. The recoil velocity is astrophysically 
significant. If it were too large, the BH can be ejected from the host galaxy with 
important consequences on BH's mass growth through hierarchical mergers. Re- 
cently, there have been a plethora of numerical [48] and analytic [49] predictions. 
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The above discussion on energy, angular momentum and linear-momentum 
can be made more rigorous and applicable to sources with non-negligible self 
gravity thanks to the work of Bondi in the 50s [50]. Let us consider a sphere S of 
volume V and radius r containing the source. Be r much larger than the source 
dimension and the gravitational wavelength (far zone). It can be proven [29,50] 
that F^, defined byl 

= / T^" d^x , (5.27) 



is a four vector with respect to Lorentz transformations. Using the relation 
dfj^r'^" = 0, we can write 



dt 



[ 9or^° d^x = ~ (f r^' n' dS . (5.28) 
Jv Js 

For yu = the above equation gives the conservation of the energy 

= -r^ j dVlT^^m, (5.29) 
for fj, — j Eq. ( 15.281 1 gives the conservation of the linear momentum 

_ = -r^ j dnr^'m. (5.30) 



6. Application to binary systems 

6.1. Inspiral waveforms at leading Newtonian order 

Let us consider a binary system with masses mi and m2, total mass M = mi + 
7712, reduced mass 1^ ~ nii 1712/(1711+7712) and symmetric mass-ratio;^ = /i/M. 
We first assume that the two bodies are rather separated and move along a circular 
orbit. In the center-of-mass frame we can write for the relative coordinates 

X{t) ^ Rcosujt, Y{t) = Rsinujt, Z{t)=Q, (6.1) 

R being the relative distance between the two bodies. The only nonzero compo- 
nents of the tensor M'^^ = fiX^ are 

Afii = i/ii?2(l + cos2tjt), (6.2) 

Af22 = ]^nR^ {I- cos2ujt) , (6.3) 
1 9 

M12 = 2^ ■ ^^'^^ 



-Note that the integration is done over a constant- time hypersurface. 
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Taking time-derivatives and plugging the above expressions in Eqs. ( 15.141 ). ( 15.151 1 
we obtain 

h+it) = r-fiR^uj^^ '- cos(2Lut) , (6.5) 

r c"^ 2 

1 4G 

h^(t) = fiR^ cos9 sm(2ujt) , (6.6) 

where we shift the origin of time to get rid of the dependence on <j> and trade the 
retarded time with t. 

For 6 = 0, i.e. along the direction perpendicular to the orbital plane, /i+ and 
hx are both different from zero, and i ih+ cx ±ie^^*'^*, thus the wave 
is circularly polarized. For 9 = 7r/2, i.e. along the orbital plane, only /i+ is 
different from zero and the wave is linearly polarized. 

The angular distribution of the radiated power is given by Eq. ( 15.241 ). It reads 

V{9) = i (1 + 6 cos^ 6* + cos^ 61) . (6.8) 

The maximum power is emitted along the direction perpendicular to the orbital 
plane, 9 = 0, where V{tt/2) = 2. Since in the case of a binary system, there 
is always a component of the source's motion perpendicular to the observation 
direction, the power radiated does not vanish in any direction. Integrating over 
the total solid angle, we obtain the total power radiated, 

^^32G^ (6.9) 
5 c5 

If we consider the binary system composed of Jupiter and the Sun, using mj = 
1.9 X 10^" g,R = 7.8 X 10^3 cm and lo = 1.68 x 10""^ Hz, we get P = 5 x 10^ 
Joules/sec. This value is tiny, especially when compared to the luminosity of 
the Sun in electromagnetic radiation Pq = 3.9 x 10^® Joules/sec. If the binary 
moves along an eccentric orbit the power radiated is [51] 

where a is the semi-major axis and e the eccentricity. Plugging in Eq. ( 16.101 ) the 
values for the Hulse-Taylor binary pulsar [52], a = 1.95 x 10^^ cm, mi = 
1.441M0, m2 = 1.383M0, e = 0.617, we get that the power radiated is 
7.35 x 10^^ Joules/sec, which is about 2% of the luminosity of the Sun in elec- 
tromagnetic radiation. 
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The emission of GWs costs energy and to compensate for the loss of energy, 
the radial separation R between the two bodies must decrease. We shall now 
derive how the orbital frequency and GW frequency change in time, using New- 
tonian dynamics and the balance equation 

^E2^^-P_ (6.11) 
dt 

At Newtonian order, i^orbit = — Gmim2/(2i?) and = GM/B?. Thus, 
R = —2/3 {Ruo) {(jj/uj'^). As long as oj/lo^ <^ 1, the radial velocity is smaller 
than the tangential velocity and the binary's motion is well approximated by an 
adiabatic sequence of quasi-circular orbits. Equation ( 16.1 11 1 implies that the or- 
bital frequency varies as 



UJ _ 96 (GMuj\ 
1> 



5/3 



2 -"l^J ' ^^-^^^ 



and the GW frequency /gw — 20", 

. 96 (GM\"^ .11/3 

where Ml ~ /i"^/^ M is the so-called chirp mass. Introducing the time to coales- 
cence T = tcoax — t, and integrating Eq. ( 16.131 ). we get 

where I.2IM0 is the chirp mass of a NS-NS binary. Equation ( 16.14b predicts 
coalescence times of ^ 17min, 2scc, Imscc, for /gw ^ 10, 100, 10'^ Hz. Using 
the above equations, it is straightforward to compute the relation between the 
radial separation and the GW frequency, we find 

Finally, a useful quantity is the number of GW cycles, defined by 

A^GW = - / uj(t)dt=- -duj. (6.16) 



Assuming ojfin 3> ti^in, we get 



A^Gw ^ 10^ ■ (6.17) 

' 1.21Moy VlOHz,' ^ ' 
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6.2. Inspiral waveform including post-Newtonian corrections 

Already in the early developments of Einstein theory, an approximation method 
called post-Newtonian (PN) method was developed by Einstein, Droste and De 
Sitter. This method allowed theorists to draw quickly several observational con- 
sequences, and within one year of the discovery of general relativity, led to the 
predictions of the relativistic advance of the perihelion of planets, the gravita- 
tional redshift and the deflection of light. 

The PN method involves an expansion around the Newtonian limit keeping 
terms of higher order in the small parameter [12, 14] 



doh 


2 






rpij 


dih 











(6.18) 



In order to be able to determine the dynamics of binary systems with a precision 
acceptable for detection, it has been necessary to compute the force determining 
the motion of the two bodies and the amplitude of the gravitational radiation with 
a precision going beyond the quadrupole formula ( 15.14b . For nonspinning BHs, 
the two-body equations of motion and the GW flux are currently known through 
3.5PN order [13]. If we restrict the discussion to circular orbits, as Eq. ( 16.12b 
shows, there exists a natural adiabatic parameter lu/uj'^ = 0[{v/c)^]. Higher- 
order PN corrections to Eq. ( 16.121 ) have been computed [13, 15], yielding (G = 
l = c) 



96 
T 



7 

k=0 



(6.19) 



where we define v^^ = {M w)^/"^ and 



WoPN 


= 1, 




1^0. 5PN 


- 0, 




WlPN 


743 

^336 


11 


^1.5PN 


= 47r + 


■ 47 Se 


W2PN 


3410J 
1814^ 


. 13 661 
i ' 2 016 



25 6m Til 

59 y 

18 
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1^X1X2 



721{e- Si)i£- S2] 



(6.20) 
(6.21) 

(6.22) 



(6.23) 



247 (Si -82)- 



(6.24) 
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Table 1 Post-Newtonian contributions to the number of GW cycles accumulated from u>i^ = n X 
10 Hz to Ljfin = u}^^'^° = 1/(6^/^ M) for binaries detectable by LIGO and VIRGO. We denote 
Ki = Si -iand^ = Si ■ S2. 



(10 + 10)Mq I (1.4 + 1.4)Mo 



Newtonian 


601 


16034 


IPN 


+59.3 


+441 


1.5PN 


-51.4+ 16. OKI XI + I6.OK2 X2 


-211 + 65.7K1 XI + 65.7K2 X2 


2PN 


+4.1 - 3.3 Ki /i2 XI X2 + l.l? XI X2 


+9.9 - 8. OKI K2X1 X2 + 2.84x1 X2 


2.5PN 


— 7.1 + 5.5 Ki XI + 5.5 K2 X2 


-11. 7 + 9. OKI XI +9.0K2 X2 


3PN 


+2.2 


+2.6 


3.5PN 


-0.8 


-0.9 



^^2.5PN = 



W3PN 



1 

672 
473 



(4159+ 15 876 1^)11 
1231 



56 

16 447 322 263 
139 708 800 



Sm Hi 
1712 



105 



451 
18"' 



4415 
4 032 



541 , 

I V 

896 

358 675 
6 048 



31811 5039 



1008 



16 



84 

(6.25) 

56198 689 



217 728 



5 605 o 



2 592 

91495 

^ 1512 



856, r_ .ji 



(6.26) 
(6.27) 



We denote £. ~ xV the Newtonian angular momentum (with X and X the 
two-body center-of-mass radial separation and relative velocity), and £ = ^/|^|; 
Si = xi m\ Si and S2 ~ X2 S2 are the spins of the two bodies (with S'1,2 
unit vectors, and < xi,2 < 1 for BHs) and 



S = Si + S2 



S = A/ 



m2 



Si_' 

mi 



(6.28) 



Finally, 6m ~ mi — 7712 and 7^; = 0.577 ... is Euler's constant. 

It is instructive to compute the relative contribution of the PN terms to the total 
number of GW cycles accumulating in the frequency band of LIGO/VIRGO. 
In Table [T] we list the figures obtained by plugging Eq. (|6.19t into Eq. (16.16b . 
As final frequency we use the innermost stable circular orbit (ISCO) of a point 
particle in Schwarzschild [/^^^° ~ 4400/(A//Afo) Hz]. 
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Fig. 4. We sketch the curvature potential as function of the tortoise coordinate r* associated to metric 
perturbations of a Schwarzschild BH. The potential peaks at the last unstable orbit for a massless 
particle (the light ring). Ingoing modes propagate toward the BH horizon, whereas outgoing modes 
propagate away from the source. 



6.3. Full waveform: inspiral, merger and ring-down 

After the two BHs merge, the system settles down to a Kerr BH and emits quasi- 
normal modes (QNMs), as originally predicted by Ref. [53,54]. This phase is 
commonly known as the ring-down (RD) phase. Since the QNMs have complex 
frequencies totally determined by the BH's mass and spin, the RD waveform is 
a superposition of damped sinusoidals. The inspiral and RD waveforms can be 
computed analytically. What about the merger? Since the nonlinearities domi- 
nate, the merger would be described at best and utterly through numerical simula- 
tions of Einstein equations. However, before NR results became available, some 
analytic approaches were proposed. In the test-mass limit, v <^1, Refs. [54,55] 
realized a long time ago that the basic physical reason underlying the presence 
of a universal merger signal was that when a test particle falls below 3A/ (the 
unstable light storage ring of Schwarzschild), the GW it generates is strongly fil- 
tered by the curvature potential barrier centered around it (see Fig. |4|i. For the 
equal-mass case v Price and Pullin [56] proposed the so-called close-limit 

approximation, which consists in switching from the two-body description to the 
one-body description (perturbed-BH) close to the light-ring location. Based on 
these observations, the effective-one-body (BOB) resummation scheme [16] pro- 
vided a first example of full waveform by (i) resumming the FN Hamiltonian, (ii) 
modeling the merger as a very short (instantaneous) phase and (iii) matching the 
end of the plunge (around the light-ring) with the RD phase (see Ref. [57] where 
similar ideas were developed also in NR). The matching was initially done using 
only the least damped QNM whose mass and spin were determined by the binary 
BH energy and angular momentum at the end of the plunge. An example of full 
waveform is given in Fig.|5] Today, with the spectacular results in NR, we are in 
the position of assessing the closeness of analytic to numerical waveforms for in- 
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Fig. 5. On the left panel we show the GW signal from an equal-mass nonspinning BH binary as 
predicted at 2.5PN order by Buonanno and Damour (2000) in Ref. [16]. The merger is assumed 
almost instantaneous and one QNM is included. On the right panel we show the GW signal from 
an equal-mass BH binai'y with a small spin Xi = X2 = 0.06 obtained in full general relativity by 
Pretorius [58] 



spiral, merger and RD. In Fig.|6l we show some first-order comparisons between 
the EOB-analytic and NR waveforms [58] (see also Ref. [59]). Similar results for 
the inspiral phase but using PN theory [13, 15] (without resummation) at 3.5PN 
order are given in Refs. [58,59]. So far, the agreement is qualitatively good, 
but more accurate simulations, starting with the BHs farther apart, are needed to 
draw robust conclusions. 

Those comparisons are suggesting that it should be possible to design purely 
analytic templates with the full numerics used to guide the patching together of 
the inspiral and RD waveforms. This is an important avenue to template construc- 
tion as eventually hundreds of thousands of waveform templates may be needed 
to extract the signal from the noise, an impossible demand for NR alone. 



6.4. Inspiral templates for data analysis 

The search for GWs from coalescing binaries with laser interferometer GW de- 
tectors is based on the matched-filtering technique, which requires accurate knowl- 
edge of the waveform (or template) of the incoming signal. As an example, in 
this section we derive the inspiral GW template in Fourier domain using the sta- 
tionary phase approximation (SPA). Those templates are currently used to search 
for inspiraling binary with LIGO/VIRGO/GEO/TAIVIA detectors. Henceforth, 
we use G = 1 = c. 
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Fig. 6. Comparison between inspiral-merger-ring-down frequency (left panel) and waveform (right 
panel) for an equal-mass BH binary with spin Xi = X2 = 0.06, as predicted at 3.5PN order in 
Ref. [16] and as obtained in a numerical simulation by Pretorius [58]. In the analytic model the 
merger is assumed almost instantaneous and three QNMs are included [58]. 'I'4 is the Weyl tensor 
proportional to the second derivative of h; we denote with tcAH the time when the apparent common 
horizon forms. 



The detector response to a GW signal is given by [24] 

hit) = h+it) F+ + h^{t) , (6.29) 

F+ = ^(l + cos^e) cos2$ cos2^'-cose sin2$ sin2*, (6.30) 
Fx = ^(1 + 008^6) cos 2$ sin 2* + cose sin 2$ cos 2*, (6.31) 

being B, $, i/' the angles defining the relative orientation of the binary with re- 
spect to the detector [24]. It is convenient to introduce the variables 



(1 + cos^ 9) 
_^4c 
4cos2eFx 



^+ [(1 + C0S2 e)^Fl+4: C0S2 e F2 ] 1/2 ' ^^'^^^ 

(6.33) 



[(1 +COs2e)2F2 +4cOs2eF2]l/2 ■ 

Noticing that Ff_ + F^ =1, we can define cos ^ = F-f and sin = Fx, and 

y^(e,$,^';e') EE [(l+cos^e)^ Fl+Acos^eF^]^/'^ , (6.34) 
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The GW signal for an inspiraling binary ( 16.291 ) can be rewritten in the simpler 
form 

h(t) = — ^(9, $, ^-^ 0) [M uj{t)]^/^ cos[2</>(t) + 200 - . (6.36) 
r 

If we are not interested in recovering the binary's orientation with respect to the 
detector (the so-called inverse problem), we can absorb f into (pQ, and average 
over the angles (9, $, 6) obtaining [24] = 16/25. 

Let us now compute the Fourier transform of the GW signal 

4-00 



Hf) = / e'-^f*h{t)dt, 



— OO 



dtA{t) 



(6.37) 



where A{t) = {2M/ R) VJ^ [M u;{t)]^/^ and (j)Gw{t) = 2(j){t) + 2$o- We 
compute the integral as follows. In Eq. ( 16.371 ) The dominant contribution comes 
from the vicinity of the stationary points in the phase. Assuming / > 0, we 
pose ip{t) = 2tt ft ~ 0GW and impose {d^p/dt)^^ = 0, that is {dcpGW /dt)^^, = 
27r / = 2-11 F{tf). Expanding the phase up to quadratic order 

^{tf)=2T,ftf~(l)Gw{tf)'T:F{tf){t-tff, (6.38) 

we get 

To compute (/>Gw(i/) and F{tf),we need to solve the following equations 

o ■ M dE , ^ 
v^ = <Pgw^, {v) - ~T{v) , (6.40) 

where E is the center-of-mass energy and the GW energy flux. A direct cal- 
culation yields 

t{v) = U + M / dv^^, (6.41) 

3 E'{v) 



(t^Gwiv) = </)c + 2 / dvw^^-^, (6.42) 



thus, 



= 2^/ - </>c - J + 2 / («,^ - i;^) --^ d^; , (6.43) 
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Using Eq. ( |6?T2] |. we have F{tf) = iIj/tt ^ (96/5) (I/tt) z^A/^/^ w"/^ and we 
obtain [60] 

= AsPA{f)e'^'^^^^\ (6.44) 

-^spa(/) = (6.45) 

TT 3 ^ 

V'spa(/) = 2Tr ftc - 0c " 7 + -r^ g- V ■0(fe/2)PN , (6.46) 

4 128 i^w? ^ 

/ A:=0 

where we denote Vf ~ [ttM f)^/^. The coefficients '0(fe/2)PN's, k = 0, . . . , 7, 
(with = 7 at 3.5PN order) in the Fourier phase are given by 

V^OPN = 1, (6.47a) 

V-csPN = 0, (6.47b) 

V'l.sPN = -167r + 4/3, (6.47d) 

/ 15293365 27145 3085 .\ „ , 

V'2PN = ^ v H - lOo- , (6.47e) 

^ V 508032 504 72 y ' ^ ' 

^2.5PN = ^f^-^z^) (l + 31og«/) , (6.47f) 
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4694215680 3 21 

15737765635 2255 tt^' 



3048192 12 



+ 



76055 2 127825 . 6848, , , 

log (4 v) , (6.47g) 

1728 1296 21 ^ ^ b' 

/ 77096675 378515 74055 



V'3.5PN = TT + ——rrr V - -T^TTT , (6.47h) 



V 254016 1512 756 



where 



^ - E(^|^ + f-)x.^., (6.48) 
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siin^ ^ sun' 

Fig. 7. We show the sky-averaged SNR for an equal-mass nonspinning binary when either the PN 
inspiral waveform (terminated at the PN ISCO) or the full NR wavefomi are included. The left panel 
uses the noise spectral density of LIGO, whereas the right panel the noise spectral density of LISA. 



In alternative theories of gravity [61], such as Brans-Dicke theory or massive 
graviton theories, the SPA phase ( 16.461 1 contains the terms 



/ altth 



if) 



128 i/w 



/ 



128 TT^DvM 



"84U.BD f 3 A2(l + z)"^ 



(6.50) 



The first term in the square brace is the contribution of dipole gravitational radi- 
ation in Brans-Dicke theory. The scalar charge of the i— th body is — acti = 
Q;(l — 2si), where ~ l/(2a;BD + 3) ^ (2ti>BD)^^ in the limit tJBD ^ 1, ands^ 
is called the sensitivity of the i— th body (a measure of the self-gravitational bind- 
ing energy per unit mass). The coefficient in the dipole term is 5 = {ai — a2)/2. 
The fact that it is dipole radiation means that it is proportional to compared 
to the quadrupole term, but the small size of S and the large current solar-system 
bound on cjbd > 40, 000 make this a small correction. The second term in the 
square brace of Eq. ( |6.50t is the effect of a massive graviton, which alters the 
arrival time of waves of a given frequency, depending on the size of the graviton 
Compton wavelength \g and on a distance quantity D, defined in Ref. [62]. It is 
possible to put interesting bounds on lo-qt) and Ag using LISA [61,62]. 

Finally, the sky-average signal-to-noise ratio (SNR) for SPA waveforms is 




1 1 



r2/3 



15 



/tin J-7/3 

'SJJ) 



-I 1/2 



(6.51) 



In the left panel of Fig. [T] we show the sky average SNRs versus total mass, for 
an equal-mass nonspinning binary at 100 Mpc using the noise-spectral density of 
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LIGO. Astrophysical observations and theoretical predictions suggest that stellar 
mass BHs have a total mass ranging between 6-3OAf0. If binary BHs of larger 
total mass exist, they could be detected by LIGO with high SNR. In the right 
panel of Fig. [T] we show the average SNR for one Michelson LISA configuration 
versus the total (redshifted) mass for an equal-mass nonspinning binary at 3 Gpc. 
The dip in the plot is due to the white-dwarf binary confusion noise. Due to the 
inclusion of merger and ring-down phases, the SNR increases considerably for 
total masses larger than 2 x IQ^Mq. 

7. Other astrophysical sources 

7.1. Pulsars 

The production of GWs from a rotating body is of considerable importance, in 
particular for application to isolated pulsars. For simplicity we shall examine the 
case the rigid body rotates around one of its principal axis. 

Let us denote by {X[ , X2, X^) the coordinates in the reference frame attached 
to the body, the so-called body frame, and introduce a fixed frame with coordi- 
nates {Xi , X2, X^), oriented such that X^ = X3. In both frames the origin of the 
axes coincides with the center-of-mass of the body. The two frames are related 
by the time-dependent rotation matrix 



where lu is the rotation frequency. The inertia tensor of a rigid body is defined by 



p being the mass density. Since any hermitian matrix can be diagonalized, there 
exists an orthogonal frame where lij is diagonal. The eigenvalues are called 
principal moments of inertia. The frame where the inertia tensor is diagonal is 
the body frame. We denote by /{^ = diag(/i, /2, /a) the inertia tensor in the 
{X[, X2, X^) coordinate system, and by lij its components in the {Xi,X2, X^) 
frame. Using the relation 7 = 7?.*/' TZ, it is straightforward to derive 
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I22 = l-i(/i-/2) cos(2t^i), (7.5) 

/33 - h, (7.6) 
/l3 = /23 = 0. (7.7) 

To get the leading-order GW signal, we need to compute the second-time deriva- 
tive of the quadrupole tensor = J d^X p(X)X' X^,i.e. Mij = -/y . We 
obtain 

A/11 = -\{h-l2)cos{2ujt), (7.8) 

M12 = -i(/i-/2) sin(2c^t), (7.9) 

M22 = i(/i -/2) cos(2c^i). (7.10) 

Thus, if the body rotates around the principal axis X3, there is a time-varying 
second mass moment only if Ii 7^/2. Plugging in Eqs. ( |5.14t , ( 15.151 ) the above 
expressions, yields 

h+ ^ /i-/2 ^ ^ ^cos 2c^0, (7.11) 

/i^ = 1 (/i - /a) cos6> sin(2cj^) . (7.12) 

r c 

The GW signal is emitted at twice the pulsar rotation frequency. It is common to 
define the ellipticity e = {h - l2)/h and = {A7r^G/c^)ll3f^^/r)e. The 
value of the ellipticity depends on the NS properties, in particular the maximum 
strain that can be supported by its crust. Pulsars are thought to form in supernova 
explosions. The outer layers of the star crystallize as the newborn pulsar cools 
by neutrino emission. Anisotropic stresses during this phase could lead to val- 
ues e < 10~^ [39], although with exotic equation of state e ~ 10^^-10^* [63]. 
Plugging numbers we find 

Using Eq. ( |5.24| i, we can compute the total power radiated 

Due to GW emission, the rotational energy of the star decreases as 
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Since the rotational energy of a star rotating around its principal axis is 



if the GW emission were the dominant mechanism for the loss of rotational en- 
ergy, the pulsar frequency should decrease as 



From electromagnetic observations one finds instead lu cx — w", where n ~ 2-3. 
Thus, the GW emission is not the major mechanism of energy loss for a rotating 
pulsar. 

If the rotation axis does not coincide with a principal axis, the pulsar motion is 
a combination of rotation around the principal axis and precession. New features 
appear in the GW signal, as discussed in detail in Ref. [64]. 

The detection of continuous, monochromatic frequency waves is achieved by 
constructing power spectrum estimators and searching for statistically significant 
peaks at fixed frequencies for very long time. If T is the observation time, the 
signal-to-noise ratio grows like \/T. The detection is complicated by the fact 
that the signal received at the detector is not perfectly monochromatic due to 
the Earth's motion. Because of Doppler shifts in frequency, the spectral lines of 
fixed frequency sources spread power into many Fourier bins about the observed 
frequency. Given the possibility that the strongest sources of continuous GWs 
may be electromagnetically invisible or previously undiscovered, an all sky, all 
frequency search for such unknown sources is very important, though computa- 
tionally very expensive [25]. 

7.2. Supenwvae 

Supernovae are triggered by the violent collapse of a stellar core which forms a 
NS or BH. The core collapse proceeds extremely fast, lasting less than a second 
and the dense fluid undergoes motions with relativistic speeds. Small deviations 
from spherical symmetry during this phase can generate GWs. 

From electromagnetic observations we know that stars with mass M > 8Mq 
end their evolution in core collapse and that 90% of them have mass between 
8-2OM0. During the collapse most of the material is ejected and if the progeni- 
tor star has a mass M < 20AIq, it leaves behind a NS. If M > 20Mq, more than 
10% of the material falls back and pushes the proto-NS above the maximum NS 
mass, leading to a BH. If the progenitor star has a mass M > AQMq, no super- 
novae takes place, the star collapses directly to a BH. Numerical simulations [65] 



(7.16) 




(7.17) 
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have predicted strains on the order of 

(7.18) 

where AE'gw = Vcff M c^, rjcs ^ 10^^. Reference [65] pointed out that GWs 
could also be produced by neutrino emission during the supernovae explosion. 
The GW signal would extend toward lower frequencies ^ 10 Hz. Moreover, the 
superposition of independent GW signals from cosmological supernovae may 
give rise to a stochastic GW background. While the estimates remain uncertain 
within several orders of magnitude, this background may become detectable by 
second-generation space-based interferometers operating around ~ 0.1 Hz [66]. 

Note that after a supernovae explosion or a collapsar a significant amount of 
the ejected material can fall back, subsequently spinning and heating the NS or 
BH. Quasi-normal modes can be excited. There is also the possibility that the col- 
lapsed material might fragment into clumps which orbit for some cycles like a a 
binary system or form bar-like structures, which also produced GW signals [67]. 



8. Cosmological sources 

In this section we want to review stochastic GW backgrounds. Depending on 
its origin, the stochastic background can be broadly divided into two classes: 
the astrophysically generated background due to the incoherent superposition of 
gravitational radiation emitted by large populations of astrophysical sources that 
cannot be resolved individually [68], and the primordial GW background gener- 
ated by processes taking place in the early stages of the Universe. A primordial 
component of such background is especially interesting, since it would carry 
unique information about the state of the primordial Universe. 

The energy and spectral content of such radiation is encoded in the spectrum, 
defined as follows 

^^GW = - /PGW(/) , (8.1) 
Pc 

where / is the frequency, pc is the critical energy density of the Universe (pc = 
3Hq/8ttG) and pew is the GWs energy density per unit frequency, i.e. 

Pgw = / dfpQY/{f). (8.2) 
Jo 

Before discussing the mechanisms that might have generated a primordial GW 
background, we review the main phenomenological bounds. 
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8.1. Phenomenological bounds 

The theory of big-bang nucleosynthesis (BBN) predicts rather successfully the 
primordial abundances of light elements. If at ieBN, the contribution of primor- 
dial GWs (or any other extra energy component) to the total energy density were 
to large, then the expansion rate of the Universe H will be too large and the 
freeze-out temperature which determines the relative abundance of neutrons and 
protons will be too high. As a consequence, neutrons will be more available and 
light elements will be overproduced spoiling BBN predictions. Detailed calcula- 
tions provide the following bound on the energy density [69] 



where N^, is the effective number of neutrino species at isBN [38]. The above 
bound extends to frequency (today) greater than ^ 10~^° Hz. More recently, 
Ref. [70] derived a similar bound by constraining the primordial GW energy den- 
sity at the time of decoupling idoc. The latter bound extends to lower frequency 
(today) ~ 10~^^ Hz, being determined by the comoving horizon size at the time 
of decoupling. 

Another important bound is the so-called COBE bound, which comes from 
the measurement of temperature fluctuations in the CMB produced by the Sachs- 
Wolfe effect. If 5T is the temperature fluctuation 



f^Gw(/) < hF 3 X 10-i«Hz < / < 10-^^Hz . (8.4) 



The lower frequency is fixed by demanding that the fluctuations should be inside 
the Hubble radius today, whereas the higher frequency by imposing that fluctu- 
ations were outside the Hubble radius at the time of the last scattering surface. 
Detailed analysis give [71] 



GWs could saturate this bound if the contribution from scalar perturbations is 
subdominant and this depends on the specific inflationary model. 

The very accurate timing of millisec pulsars can constrain Oqw. If a GW 
passes between us and the pulsar, the time of arrival of the pulse is Doppler 
shifted. Many years of observation yield to the bound [72] (see also Ref. [10]) 




(8.3) 





(8.5) 
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(8.6) 
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8.2. Gravitational waves produced by causal mechanisms 

Two features determine the typical frequency of GWs of cosmological origin 
produced by causal mechanism: (i) the dynamics, which is model dependent and 
(ii) the kinematics, i.e. the redshift from the production time. Let us assume that 
a graviton is produced with frequency at time during matter or radiation 
era. What is the frequency today? We have 

/ = /*-, (8.7) 
ao 

assuming that the Universe evolved adiabatically, so gs{T^:) = gs{To) Tq Oq, 

using To — 2.73 K and gs{To) = 3.91, where gs{Tt) is the number of degrees 
of freedom at temperature , we get 

What is /,? Since the size of the Hubble radius is the scale beyond which causal 
microphysics cannot operate, we can say that from causality considerations, the 
characteristic wavelength of gravitons produced at time t» is H^^ or smaller. 
Thus, we set 

If the GW signal is produced during the radiation era, 

H':-^P^8^'9*T^^-^, (8.10) 



90 Af2j 



and we find 



f^lO-^U^] (Ji^Y^'R.. (8.11) 

From the above equation, we obtain that millisec pulsars can probe physics at the 
^ MeV scale, LISA at the TeV scale, ground-based detectors at 10^-10'' TeV 
and detectors in the GHz bandwidth [73] at GUT or Planck scale. 

As an application of GWs produced by causal mechanisms, let us consider 
GW signals from first-order phase transitions. In the history of the Universe 
several phase transitions could have happened. The quantum-chromodynamics 
(QCD) phase transition takes place at T 150 MeV. Around T ^ 100 GeV 
the electroweak (EW) phase transition happens (and SU{2) x U{1) breaks to 
C/em(l) through the Higgs mechanism). Let us assume that V{(j), T) is the poten- 
tial associated to the phase transition, where is the order (field) parameter. As 
the Universe cools down, the true false vacuum are separated by a potential 
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barrier. Through quantum tunneling, bubbles of true vacuum can nucleate. The 
difference of energy between the true and false vacua is converted in kinetic en- 
ergy (speed) of the bubble's wall. In order to start expanding, the bubbles must 
have the right size, so that the volume energy overcomes the shrinking effects of 
the surface tension. The relevant parameter is the nucleation rate T = Tq e^*. If 
r is large enough, the bubbles can collide within the Hubble radius, and being 
the collision nonsymmetric, they produce GWs [74]. The parameter (3 fixes the 
frequency at the time of production. 

So, we can write for the parameter e defined in Eq. ( 18.91 ). e — H^,/ f^, — H^,/(3 
and the peak of the GW spectrum occurs at [75, 77, 79] 

'(^)(Tiv)(li)"°-- 

where T, is the temperature at which the probability that a bubble is nucleated 
within the horizon size is 0(1). 

In the case of EW phase transitions, we have ~ 10^-10^ and T, ~ 10^ 

GeV, thus /peak — 10^^-10^'^ Hz, which is in the frequency band of LISA. The 
GW spectrum can be computed semianalytically, it reads [75, 77, 79] 

where a is the ratio between the false vacuum energy density and the energy 
density of the radiation at the transition temperature T*; k quantifies the fraction 
of latent heat that is transformed into bubble-wall kinetic energy and Vf, is the 
bubble expansion velocity. 

Nonperturbative calculations obtained using lattice field theory have shown 
that there is no first-order phase transition in the standard model if the Higgs 
mass is larger than Mw- In minimal supersymmetric standard model, if Higgs 
mass is ^ 110-115 GeV, the transition is first-order but /iq ^GW ^ 10^^^ at 
/ ~ 10 mHz. In next-to-minimal supersymmetric standard models, there ex- 
ist regions of the parameter space in which ho ficw ~ 10^^^-10~^° at / ~ 10 
mHz. Note that for frequencies 10^^-3 x 10^"^ Hz the stochastic GW background 
from white-dwarf binaries could hide part of the GW spectrum from first-order 
phase transitions. More recently, Ref. [78] pointed out that new models of EW 
symmetry breaking that have been proposed have typically a Higgs potential dif- 
ferent from the one in the Standard Model. Those potentials could lead to a 
stronger first-order phase transition, thus to a more promising GW signature in 
the milliHz frequency range. 

A stochastic GW background could be also produced during a first-order 
phase transition from turbulent (anisotropic) eddies generated in the background 
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fluid during the fast expansion and collision of the true-vacuum bubbles [76, 77, 
79]. In the next-to-leading supersymmetric standard model there exist regions 
of the parameter space where [77] /iq^Iqw ~ 10~^° with peak frequency in 
the mHz. Reference [80] evaluated the stochastic GW background generated by 
cosmic turbulence before neutrino decoupling, i.e. much later than EW phase 
transition, and at the end of a first-order phase transition if magnetic fields also 
affect the turbulent energy spectrum. The observational perspectives of those 
scenarios are promising for LISA. 

Long time ago Turner and Wilczek [81] pointed out that if inflation ends with 
bubble collisions, as in extended inflation, the GW spectrum produced has a peak 
in the frequency range of ground-based detectors. Subsequent analyses have 
shown that in two-field inflationary models where a field performs the first-order 
phase transition and a second field provides the inflationary slow rolling (so- 
called first-order or false vacuum inflation [82]), if the phase transition occurs 
well before the end of inflation, a GW spectrum peaked around 10-10"^ Hz, can 
be produced [83], with an amplitude large enough, depending on the number of 
e-foldings left after the phase transition, to be detectable by ground-based inter- 
ferometers. A successful detection of such a spectrum will allow to distinguish 
between inflation and other cosmological phase transitions, like QCD or EW, 
which have a different peak frequency. 

Another mechanism that could have produced GWs is parametric amplifica- 
tion after preheating [84]. During this phase classical fluctuations produced by 
the oscillations of the inflaton field (f> can interact back, via parametric resonance, 
on the oscillating background producing GWs. In the model where the inflaton 
potential contains also the interaction term ^ 0^ X^, X being a scalar field, the 
authors of Ref. [84] estimated f^GW ^ 10~^^ at /min ~ 10^ Hz, while in pure 
chaotic inflation Qgw < 10"" at /min ^ 10"* Hz. [See Fig. 3 in Ref. [84] for 
the GW spectrum in the range 10^-10^ Hz and Refs. [85] for a recent reanalysis.] 
Unfortunately, the predictions lay in the frequency range where no GW detectors 
exist, but some proposals have been made [73] 

8.3. Gravitational waves produced by cosmic and fundamental strings 

Topological defects could have been produced during symmetry-breaking phase 
transitions in the early Universe. Since the 80s they received significant attention 
as possible candidates for seeding structure formation. Recently, more accurate 
observations of CMB inhomogeneities on smaller angular scales and compatibil- 
ity with the density fluctuation spectrum on scales of 100 h^^ Mpc, restrict the 
contribution of topological defects to less than 10%. 

Cosmic strings are characterized by a single dimensional scale, the mass-per- 
unit length fi. The string length is defined as the energy of the loop divided by 
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fi. Cosmic string are stable against all types of decay, except from the emission 
of GWs. Let us assume that a network of cosmic strings did form during the 
evolution of the Universe. In this network the only relevant scale is the Hubble 
length. Small loops (smaller than Hubble radius length) oscillate, emit GWs and 
disappear, but they are all the time replaced by small loops broken off very long 
loops (longer than Hubble radius). The wavelength of the GW is determined by 
the length of the loop, and since in the network there are loops of all sizes, the GW 
spectrum is (almost) flat in a large frequency band, extending from / ~ 10~® Hz 
to / - 10^° Hz. 

In 2001, Damour and Vilenkin [86] (see also Ref. [87]), worked out that strong 
bursts of GWs could be produced at cusps (where the string reaches a speed very 
close to the light speed) and at kinks along the string loop. As a consequence 
of these bursts the GW background emitted by a string network is strongly non- 
Gaussian. The most interesting feature of these GW bursts is that they could be 
detectable for a large range of values of Gfi, larger than the usually considered 
search for the Gaussian spectrum. GW bursts can be also produced by fundamen- 
tal strings, as pointed out in Ref. [88]. For a detailed analysis of the prospects 
of detecting the stochastic GW background and the GW bursts with ground and 
space-based detectors, and msec pulsars see Refs. [10, 89]. 

A GW burst emitted at the cusp of cosmic or fundamental strings can be de- 
tected using matched filtering. In Fourier domain the signal is [89] 



where A ~ {GfiL^^'^)/r, /e is determined by the size L of the feature that pro- 
duces the cusp, but also by the low-frequency cutoff frequency of the detector 
and fh ~ 2/{9'^L), 9 being the angle between the fine of sight and the direction 
of the cusp. 

8.4. Gravitational waves produced during inflation 

The amplification of quantum vacuum fluctuations is a common mechanism in 
quantum field theory in curved space time [90] . In the 70s Grishchuk and Starobin- 
sky [91] applied it to cosmology, predicting a stochastic GW background which 
today would span a very large frequency band 10~^^-10^" Hz. Henceforth, we 
shall compute the GW spectrum using semiclassical arguments, and refer the 
reader to Refs. [92, 93] for more detailed computations. 
The background field dynamics is described by the action 



M/)-^l/r'/'e(/,-/)e(/-/e), 



(8.14) 




(8.15) 
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We assume an isotropic and spatially homogeneous Friedmann- Lamaitre-Robertson- 
Walker metric with scale factor a, 

ds^ = -dt^ + a'^ dx^ ^ gf.r.dx'' dx" . (8.16) 

For what we have learned in previous lectures, we can derive the free-linearized 
wave equations for the TT metric perturbations Jg^j, — h^j^^, with /i^o = 0, /i^ = 
0, hi^.^ = and (5T^ = 0, obtaining 

^^d^{^g''''d,)hl =0, (8.17) 

where we disregard any anisotropic stresses. Introducing the conformal time 77, 
with di] — dt/a{t), we can write 

hi{v,^) = V8^ ^ 5^/ii^(r;)e^'^-e,f^"(n), (8.18) 

A=+,x k 

ef-'{n) being the polarization tensors. Since we assume isotropic and spatially 
homogeneous metric perturbations /ik = hk, and each polarization mode satisfies 
the equation 

h'^iij) + 2- h'kirj) + hkiv) = , (8.19) 
a 

where we denote with a prime the derivative with respect to the conformal time. 
By introducing the canonical variable ipk {r}) = o-hk [i]) W6 can recast Eq. ( |8.19l l 
in the much simpler form 

i^kiv) + [fc' - Uiv)] Vfc = , U{rj) - — , (8.20) 

a 

which is the equation of an harmonic oscillator in a time-depedent potential U{ri). 
We want to solve the above equation and study the properties of the solutions. 

For simplicity, we consider a De Sitter inflationary era, a = —1/ {Hus v) 
a" /a ~ '^/rf', and make the crude assumption that the De Sitter era is followed 
instantaneously by the radiation era, 0(77) = (77 — 277*)/(77ds '?*) anda"(77) = 0. 

If A:^ ^ U{ni), i.e. kij S> 1 or a//c ^ jaTyj = iJ^g, the mode is inside the 
Hubble radius or (in jargon) is over the potential barrier U{rj) and the (positive 
frequency) solution reads 
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thus hk decreases while inside the Hubble radius. If fc^ <C U{ri), i.e. krj ^ 1, 
a/k ^ \ari\ = H^g, the mode is outside the Hubble radius or (in jargon) under 
the potential barrier. In this case the solution reads 



Ak + Bk 



drj 



Ak+B,lJ^y (8.22) 



Since during the De Sitter era the scale factor gets larger and larger, the term 
proportional to B in the above equation becomes more and more negligible. Thus, 
the perturbation hk remains (almost) constant while outside the Hubble radius. 
So, the longer the tensorial-perturbation mode remains outside the Hubble radius, 
the more it gets amplified (with respect to the case it stayed always inside the 
Hubble radius). During the RD era, the mode is again under the barrier, and the 
solution is 

^fc = ak e-'^" + Pk 6+"^" , (8.23) 

and contains both positive and negative modes. Even normalizing the initial state 
to positive frequency mode, the final state is a mixture of positive and frequency 
modes. In a quantum field theory language, such a mixing represents a process 
of pair production from vacuum. The coefficients ak and /3k are called Bogoli- 
ubov coefficients and can be obtained imposing the continuity of the tensorial 
perturbation and its first time derivative at the transition between cosmological 
phases. 

The intensity of the stochastic GW background can be expressed in terms of 
the number of gravitons per cell of the phase space Uf with / = |k|/ (27r). For 
an isotropic stochastic GW background pcw = 2 J cfik/ {2tt)^ [k Uk), thus 

f^Gw(/) = -16^'n//4. (8.24) 

Pc 

where n/ = 

The stochastic GW spectrum produced during slow-roll inflation, decreases 
as in the frequency window 10^^®-10^^^ Hz, and then slowly decreases 
up to a frequency corresponding to modes whose physical frequency becomes 
less than the maximum causal distance during the reheating phase (which is of 
order of a few GHz). Its magnitude depends on both the value of the Hubble 
parameter during inflation and a number of features characterizing the Universe 
evolution after the inflationary era — for example tensor anisotropic stress due 
to free-streaming relativistic particles, equations of state, etc. [93,94]. An upper 
bound on the spectrum can be obtained from the measurement of the quadrupole 
anisotropy of the CMB, as seen in Sec. 18.11 Since for a generic slow-roll infla- 
tionary model the spectrum is (weakly) decreasing with frequency this implies an 
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upper bound ficw ^^-^ 5 x 10~^^ at frequencies around / ^ 100 Hz , where 
ground-based detectors reach the best sensitivity. The spectrum predicted by the 
class of single-field inflationary models is then too low to be observed by ground- 
based and also space-based detectors. It is therefore evident that a background 
satisfying the bound imposed by the observed amount of CMB anisotropics at 
large scales could be detected by GW interferometers provided that its spectrum 
grows significantly with frequency. This could happen in bouncing Universe 
cosmologies, such as pre-big-bang scenario [95,96], the ekpyrotic models [97] 
(although the amplitude of the GW spectrum is too low to be observable) and in 
quintessential inflation [98]. 

Finally, a stochastic GW background can be detected by correlating two GW 
interferometers [99]. The upper limit on a flat-spectrum set by the LIGO Sci- 
entific Collaboration is /iqSIgvv — 6.5 x 10^^ in the frequency band 70-156 
Hz [28]. For frequency-independent spectra, the expected upper limit for the 
current LIGO configuration is /i§ ilcvv < 5 x 10^^, while for advanced LIGO 
project is flow ^ 8 x 10^^. 
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